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The spatial distribution of partons in the proton can be probed in suitable exclusive scat- 
tering processes. I report on recent performance estimates for parton imaging at a proposed 
Electron-Ion Collider. 

Deeply virtual Compton scattering (DVCS) and exclusive meson production in lepton-proton 
collisions offer unique possibilities for determining the spatial distribution of quarks, antiquarks 
and gluons as a function of their longitudinal momentum inside the proton. Such "tomographic 
images" of the proton can provide insight into key aspects of QCD dynamics, such as the 
interplay between sea quarks and gluons, the relation between sea and valence quarks, and the 
orbital angular momentum carried by these partons. 

To obtain such images one needs to measure the transverse momentum transfer A to the 
proton. A Fourier transform then gives the distribution in impact parameter b, which is the 
position of the struck parton in the plane transverse to the proton direction of movement. In 
the scattering amplitude the longitudinal momentum fraction x of the parton is integrated over, 
with typical values around \xb for DVCS and \xy = \xb (1 + M 2 /Q 2 ) for the production of 
a meson with mass M. The reconstruction of a joint density in x and b can be envisaged in the 
framework of generalized parton distributions (GPDs) by making use of their evolution in the 
resolution scale, which is given by Q 2 in DVCS and by Q 2 + M 2 in meson production. This 
requires precise data in a very wide range of xb and Q 2 . 

Measurements from HERA suggest that at x around 10~ 3 gluons have a more narrow impact 
parameter distribution than sea quarks. Information about the spatial distribution of valence 
quarks can be inferred from electromagnetic form factors and from lattice calculations, and 
a direct study of quarks with large momentum fraction x will be possible with the 12 GeV 
upgrade at Jefferson Lab. The planned DVCS measurement by COMPASS will give us a first 
glimpse into the x region between 10 _1 and 10~ 2 . To realize the full physics potential of parton 
imaging will, however, require a new facility. Here I report on a study of parton imaging at a 
proposed electron-ion collider (EIC) [T]. For complementary information see [2J[3]. 

Pseudo-data for ep — > epj have been generated according to a GPD model that reproduces 
the existing DVCS measurements of HI and ZEUS. Technically, the unpolarized distributions H 
for gluons and for sea quarks are modeled with two SO (3) partial waves each, as described in [3] 
and in Sec. 3 of [I]. The t dependence of the distributions is exp[t (-j +o/ log ^)] , with different 
slopes B for gluons and sea quarks and with a small shrinkage parameter a' as suggested by 
HERA measurements. In DVCS, the invariant momentum transfer t to the proton and its 
transverse component are related by —t = (A 2 + x 2 B m 2 p )/{l — x B ), where m p is the proton 
mass. In the simulation [2] acceptance cuts for the final-state electron, photon and proton have 
been imposed, and the data have been smeared for the expected resolution in xb, Q 2 and t. 
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Figure 1: Left: simulated DVCS cross section at EIC for a bin in xb and Q 2 . Right: average 
squared impact parameter obtained from der/d^DVCS for different bins in xb and Q 2 . 
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An assumed 5% uncorrelated systematic error has been added in quadrature to the statistical 
error in the cross section. The cross section da/dt for DVCS (j*p — ► jp) is obtained after 
subtraction of the cross section for the Bethe-Heitler process, with an uncertainty of 3% taken 
on the latter. Beam energies are E e — 20 GeV and E p = 250 GeV. Statistical errors are for 
an integrated luminosity of 10 fb _1 for the |t| range from 0.03 GeV 2 to 1 GeV 2 and for 100 fb _ 
for \t\ above 1 GeV 2 . Here and in the following the uncertainty in the overall luminosity is not 
included in the errors, because it does not affect the form of the spectra which are at the center 
of our interest. 

The left panel of Fig. [T] shows the result- 
ing t spectrum for DVCS in a bin of xb and 
Q 2 , together with an exponential fit. From the 
t slope one obtains the average squared im- 
pact parameter (b 2 ) for the particular combi- 
nation of quarks, antiquarks and gluons "seen" 
in DVCS. The right panel of Fig. Q] shows that 
with the expected accuracy one can resolve the 
separate dependence of (b 2 ) on Q 2 and xb, 
which has never been possible so far. In the 
model used for generating the data, both de- 
pendences are small logarithmic effects, which 
reflect both perturbative and non-perturbative 
dynamics of sea quarks and gluons in the pro- 
ton. From da/dt one obtains the Compton 
scattering amplitude |_4 7 * P ^ 7P |, whose Fourier 
transform w.r.t. A gives an impact parameter 
distribution F{xb, b) of partons with momen- 
tum fraction of order \xb- Figure [2] shows 
that precise imaging is possible in a wide range 
of b. Of particular interest at large b is the prediction of an exponential falloff with an x de- 
pendent slope of typical size l/(2m ff ) « 0.7 fm [S]. Due to partons inside virtual pions, which 
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Figure 2: Impact parameter distribution in 
DVCS obtained from the left panel of Fig. [TJ 
The band reflects the uncertainty in fitting 
da/dt and in extrapolating it to the unmea- 
sured t region as specified in Sec. 3.6.1 of pQ. 
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because of their small mass can fluctuate to large distances, this predicted behavior is a conse- 
quence of chiral symmetry breaking in QCD and awaits experimental verification. 




Figure 3: Impact parameter densities of unpolarized sea quarks in an unpolarized proton (left) 
and in a proton polarized along the x axis (right), obtained from a fit to pseudo-data for the 
DVCS cross section and the transverse proton spin asymmetry [3J [3] . 



With transverse proton polarization one 
gains access to distributions E, which carry 
characteristic information about orbital an- 
gular momentum of partons in the proton. 
Transformed to b space, these distributions de- 
scribe how the impact parameter distribution 
of partons is shifted sideways in a transversely 
polarized proton. Pseudo-data for the DVCS 
cross section and for the transverse proton spin 
asymmetry have been generated [2] assuming a 
model of E for sea quarks and for gluons of the 
same type as the model of H described above. 
Parameter are chosen to satisfy the positiv- 
ity requirements on E. The simulation is for 
E e = 20 GeV and E p = 250 GeV, with errors 
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Figure 4: As Fig. [31 but for unpolarized gluons 
in an unpolarized proton. 



for an integrated luminosity of 100 fb _ and 
for 80% transverse proton polarization mea- 
sured with 5% accuracy. Systematic uncertainties of 5% are added in quadrature. As shown 
in Fig. [31 a fit of the generated data for the DVCS cross section and the transverse proton 
spin asymmetry allows the simultaneous extraction of the GPDs H and E for sea quarks. The 
fitted data covers the kinematic range 3.2 GeV 2 < Q 2 < 17.8 GeV 2 , 10" 4 < x B < 10~ 2 and 
0.03 GeV 2 < \t\ < 1.5 GeV 2 . Thanks to the effect of logarithmic scaling violations, even the 
extraction of H for gluons is possible in this fit, as shown in Fig. 21 whereas the errors on E for 
gluons are very large. 

Selective information about gluons can be obtained in exclusive J/ VP production. We have 



generated pseudo-data for e p 



5 (e + e )p with a version of PYTHIA modified 



to describe the HI and ZEUS measurements of this process. We assume an acceptance in 
pseudorapidity of rj < 5 for all final-state leptons and an acceptance for the recoil proton as in 
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the DVCS studies above. The t spectra and corresponding impact parameter distributions in 
Fig. [5] show that EIC can provide accurate images of gluons in the proton over two orders of 
magnitude in x. 



0.0016 <x v < 0.0025 15.8 GeV 2 < Q 2 + < 25.1 GeV 2 0.16 <x v < 0.25 10.5 GeV 2 < Q 2 + M 2 ,^ < 15.8 GeV 2 
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Figure 5: Simulated cross sections for 7*p — > J/^t p and the corresponding impact parameter 
distributions for gluons. The left panels are for E e = 20 GeV, E p = 250 GeV and the right 
panels for E e = 5 GeV, E p = 100 GeV, with 10 fb — integrated luminosity in both cases. 
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